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- - - By LDA+U method with spin-orbit couphng (LDA+U+SO) the magnetic state and electronic 

structure have been investigated for plutonium in 5 and ct phases and for Pu compounds: PuN, 
■ PuCoGas, PuRh2, PuSi2, PuTe, and PuSb. For metallic plutonium in both phases in agreement 

with experiment a nonmagnetic ground state was found with Pu ions in f^ configuration with zero 
values of spin, orbital, and total moments. This result is determined by a strong spin-orbit coupling 
. in 5/ shell that gives in LDA calculation a pronounced splitting of 5/ states on f"'^'^ and /^''^ 

' subbands. A Fermi level is in a pseudogap between them, so that /^^^ subshell is already nearly 

completely filled with six electrons before Coulomb correlation effects were taken into account. 
The competition between spin-orbit coupling and exchange (Hund) interaction (favoring magnetic 
^Sj ' ground state) in 5/ shell is so delicately balanced, that a small increase (less than 15%) of exchange 

interaction parameter value from Jh ~ 0.48 eV obtained in constrain LDA calculation would result 
in a magnetic ground state with nonzero spin and orbital moment values. For Pu compounds 
^ ■ investigated in the present work, predominantly /® configuration with nonzero magnetic moments 

was found in PuCoGas, PuSi2, and PuTe, while PuN, PuRh2, and PuSb have configuration 
with sizeable magnetic moment values. Whereas pure jj coupling scheme was found to be valid for 
^ . metallic plutonium, intermediate coupling scheme is needed to describe 5/ shell in Pu compounds. 

The results of our calculations show that exchange interaction term in the Hamiltonian must be 
treated in a general matrix form for Pu and its compounds. 
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PACS numbers: 71.27.-|-a, 71.70.-d, 71.20.-b 



I. INTRODUCTION 



' In actinide elements Coulomb interaction of 5/ electrons is of the same order as a band width and spin-orbit coupling. 
^ , Their interplay results in a complicated physics of 5/ compoundsi where both theorists and experimentalists still have 
' attractive plenty of work. Among other actinides, plutonium element seams to be most intriguing f^i^ despite numerous 
works and 60 years passed since the Pu discovery^*^ 

Electronic properties of plutonium show an exceptional example of the system with 5/ electrons on the edge between 
\ localization and itinerancy^*^ Experimental worki gave jj coupling for Pu 5/ electrons or an intermediate coupling 

■ close to jj typei^ Theoretical model for intermediate coupling scheme was suggested for PuSb by Cooper et al^ 

■ Many experiments have been carried out to shed light on plutonium electronic structure. For (5-plutonium, most in- 
"j--, vestigated of all metallic Pu phases, photoemission spectroscopy (PES) revealed special features of valence bandiS*ii*i^ 

in particular a sharp peak at the Fermi level. X-ray, high- resolution ultraviolet,^^ and resonant photoemissioni^ spec- 
^ troscopy measurements accompanied with Pu 4/ core-level spectra show more localized character of 5/ electrons 
I in (5-Pu phase comparing with a phase. Sharp peak at the Fermi level is indicative of strong many-particle nature 
■ O I and heavy-fermion behavior of 5/ electrons. The electron mass enhancement is usually characterized by the inferred 
^ ■ Sommerfeld coefficient 7 = 50 mJ mol~^ for (5-Pui^ versus 7 = 17 mJ mol~^ for a-Pum^ Phonon dispersions 
and elastic moduli of 6-Pu showed a number of anomalies including softening of the transverse [111] modesji^ 
suggesting unusual electronic structure and phonon-electron interaction for 5/ electrons in plutonium (modern calcu- 
lation results^ qualitatively predicted experimental data). No evidence for ordered or disordered magnetic moments 
in plutonium in both 6 and a phases was found in experiments iSii^i 

Ab-initio calculations of the electronic structure for plutonium have been done by many authors. Standard local 
, density approximation (LDA)S as well as generalized gradient approximation (GGA)^'^ have been implemented. Full- 
potential linearized augmented plain- wave method (FLAPW) with GGA-correction and Gaussian-type orbitals-fitting 
function (LCGTO-FF) methods were used to reproduce zero-pressure experimental volumes and bulk module of Pu 
in the row of all actinides. The obtained results show better agreement with experiment than, for example, LMTO- 
GGA method^^ but fails handling with Pu. Even with proper second variational approach to 6p states^Si experimental 
volume of S-Pn is still unattainable!^ Another limitation of all these methods, as it was realized later, is inability to 
reproduce correctly electronic structure of plutonium due to neglect of Coulomb correlation effects. 

To improve LDA approximation LDA+U mcthod^^. was used, which explicitly takes into account strong Coulomb 
repulsion between 5/ electrons by adding Hubbard-like term to the LDA Hamiltonian. The calculation of electronic 
structure for fcc-Pu by LDA-f U method with GGA correction was performed in Ref. Authors showed that for 
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the standard set of Slater integrals^ii^ and Coulomb parameter U = 4.0 eV significant improvement was achieved for 
calculated value of equilibrium volume comparing with experimental data. The resulting band structure suggests that 
Pu has five localized / electrons with substantial orbital and spin moments values. Another LSDA+U calculation 
with GGA correction gave similar results for different Coulomb parameter values iS 

In contrast to experimental data, clearly showing nonmagnetic state of plutonium;2ii*2i antiferromagnetic ground 
state was found for fcc-Pu in many calculations i^^i^^iSSi^l The fully relativistic spin-polarized calculation for fcc- 
plutonium phase was performed in Ref . |3^. In another work (Ref. 'SO*) it was shown that 'switching off' the hybridiza- 
tion makes better agreement with experiment for equilibrium volume in (5-Pu. For augmented plane waves basis 
(FLAPW),**- Kutepov et alm^ using the fully relativistic spin-polarized (RSP) method with GGA approximation ob- 
tained large orbital and spin moments. Antiferromagnetic ordered state was found lower in total energy of the system 
comparing with nonmagnetic one for both a- and 5-Pu phases4iii^ 

Another scheme was proposed by Eriksson et alm^ Plutonium 5/ electrons were divided into localized and delocalized 
ones. Self-interaction correction SIC-LSDA method^^i^^i^Sdi was used to calculate volumes and total energies of Pu 
ion in different configurations (with different valencies^^^ of Pu ion). In contrast to experimentSfi L5-coupled 
ground state was found to be lower in the total energy for any configuration of Pu ion^^ and configuration 
energetically preferable. For many binary Pu compounds SIC-LSDA calculations showed that experimental data are 
better described by the model of coexisting localized and delocalized electrons.''^ In particular, valency -1-5 was found 
for Pu ion in PuN and -1-3 in PuTe and PuSb4^ Pu ion in plutonium dioxide by the same approach was obtained in 
/■* configuration42i^ Similar hybrid schemes, disordered local moment (DLM) method^ ^ gave qualitative agreement 
with the photoemission spectrum for S-Pum^ 

Experiments suggest the presence of significant correlation effects for 5/ electrons in plutonium. Dynamical Mean- 
Field TheorySi^^iS (DMFT) is a powerful tool in studying such effects. Merging of LDA-based methods with DMFT 
gave a new calculation scheme - LDA-I-DMFT methodi^SiSiS Its application to plutonium problem^®'^- gave a better 
agreement with photoemission experiments and a double minimum curve of the total energy as a function of volume. 

In this paper we present the results of calculation by LDA-I-U with spin-orbit coupling (LDA-I-U-I-SO) method of 
electronic structure and magnetic properties for metallic plutonium in a and S phases and series of Pu compounds, 
namely PuN, PuCoGas, PuRh2, PuSi2, PuTe, and PuSb. Exchange interaction (spin polarization) term in the 
Hamiltonian was implemented in a general nondiagonal matrix form. This form is necessary for the correct description 
of 5/ electrons for the case of jj coupling and intermediate coupling schemes which we have found to be valid in pure 
metallic Pu and its compounds. For metallic Pu in both S and a phases, our calculations gave nonmagnetic ground 
state with six / electrons in fully occupied j — 5/2 sub-shell and spin, orbital, and total moments values equal to zero. 
This result is in agreement with experimental data for metallic plutoniumjSi but we have found magnetic ground 
state for all Pu-based compounds investigated in this work. For the present investigation we have chosen plutonium 
compounds with formal Pu ion valency -1-3 (PuN and PuSb) and +2 (PuTe) and also intermetallic compounds where 
valence state of Pu ion is not obviously defined. All calculated Pu-compounds possess magnetic moments on Pu 
ions with or mixed configurations. We were interested mainly in magnetic state problem and discussion of 

equilibrium volume and volume transition in plutonium and its compounds could be found elsewherei^ A strong 
competition between spin-orbit coupling and exchange interactions with a delicate balance was found to determine 
magnetic state of Pu and its compounds. 

The organization of this paper is as follows. Sec. describes LDA+U with spin-orbit coupling (LDA-(-U+SO) 
method used in the present work. In Sec. IIIII results of calculations for metallic 6- and a-Pu phases are discussed in 
detail. Then Sec. IIVI describes applications of LDA-I-U-I-SO method to various plutonium compounds, including PuN 
in Subsec. Hvn PuCo Gas in Subsec. UVBl PuRha in Subsec HvOl PuSia in Subsec UVDl PuTe in Subsec. HVEl and 
PuSb in Subsec. QV Sec. IVI summarizes the paper. 

II. METHOD 

A strong spin-orbit coupling in actinides together with magnetism and Coulomb interactions rise a problem which 
was not present for materials without 5/ electrons. For 3d and 4/ elements compounds relatively weak spin-orbit 
coupling results in Russell-Saunders coupling [LS coupling) scheme with S and L operators well defined. Then the 
basis of LS orbitals, which are eigenfunctions of both spin S and orbital moment L operators, is a good choice. In this 
case it is possible to define quantization axis in the direction of spin moment vector so that occupation and potential 
matrices will be diagonal in spin variables. When spin-orbit coupling is stronger than exchange (Hund) interaction, jj 
coupling scheme is valid with a well defined total moment J, but not spin S and orbital L moments. In this case, the 
basis of eigenfunctions of total moment operator {jrrij} is a best choice. The matrix of spin-orbit coupling operator 
is diagonal in this basis but not the exchange interaction (spin-polarization) term in the Hamiltonian. 

The situation for 5/ electrons is more complicated. The strengths of spin-orbit coupling and exchange interaction are 
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comparable so that both aforementioned hmits: LS coupHng and jj coupling schemes are not valid and intermediate 
coupling scheme is needed to describe 5/ shell in actinides. In this case occupation matrix is diagonal neither in 
{LS*}, nor in {jnij} orbital basis and both terms in the Hamiltonian: spin-orbit coupling and exchange interaction, 
must be taken in a general nondiagonal matrix form. 

In LDA+U method^^ the energy functional depends, in addition to the charge density p(r), on the occupation 
matrix n^^, (LDA+U method in general nondiagonal in spin variables form was defined in Ref. l6l|) : 

E^^^+^[p{r), {n}] = E^^^ipir)] + E"[{n}] - E,4{n}] (1) 

where p{r) is the charge density and E^-'^^[p{r)] is the standard LDA functional. The occupation matrix is defined 
as: 



2 /"^f 



/mG-„,(£;)di? (2) 



where G^^^, (E) — {inlms \ [E — Hj^da^jj) ^ \ inlm's') are the elements of the Green function matrix in local orbital 
basis set (i - atomic index, n - principal quantum number, Im ~ orbital quantum numbers, and s — spin index). In 
the present work this basis set was formed by LMT-orbitals from LMTO methodj^ In Eq. Q Coulomb interaction 
energy _E'^[{ri}] term is a function of occupancy matrix nt° 
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E^'iin}] = sAi^^^" I Vee I m',m"')< 

(3) 

-(to, to" I Vee I m"\m'))nf^'„^,nf^J,„^,„} 

where Vee is the screened Coulomb interaction between the nl electrons. Finally, the last term in Eq. ^ correcting 
for double counting is a function of the total number of electrons in the spirit of LDA which is a functional of total 
charge density: 

Edc[{n}] = ^UN{N - 1) - ^JN{N - 2) (4) 

where N — Tr{n^^,) is a total number of electrons in the nl shell. U and J are screened Coulomb and Hund exchange 
parameters which could be determined in constrain LDA calculations The latter (J) in this paper is everywhere 
denoted as Jh not to be confused with total moment operator J. The screened Coulomb interaction matrix elements 
{m,m" I Vee I rn',m"') could be expressed via parameters U and Jh (see Ref. 11^). 

The functional Eq. Q defines the effective single-particle Hamiltonian with an orbital dependent potential added 
to the usual LDA potential: 

Hlda+u^Hlda+ Y1 \inlms)V^'^,{inlm's'\ (5) 

ms^ni' s' 

Vmm' = 5,s' Y.m",m'"{{^^^" I ^ee I m' , iv! " ) f^^'^ , + 
{{m,m" I Vee I to', to'") - (to, to" I Vee \ to'" , to') }- 

(6) 

(1 - I Vee I TO'",TO'))n^?,^,„ 

-U{N-^) + ij(7V- 1). 

In this paper we used method abbreviated as LDA-f U-I-SO which comprises non-diagonal in spin variables LDA+U 
Hamiltonian Eq. © with spin-orbit (SO) coupling realized on the base of TB-LMTO-ASA codet^Si 

Hlda+u+so = Hlda+u + Hso, 
Hso = A • L • S. 
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In LS basis spin-orbit (SO) coupling matrix has diagonal {Hso)m' m ^ '^^^^ off-diagonal {Hso)]^' m ^'^'^ i^so)m' m 
non-zero matrix elements (in complex spherical harmonics); 

{Hso)ln',m - ^Vil+^)il-^+^)iSm'-l,.n) (8) 

where s = +1/2, -1/2^ 

The exchange interaction is originally a part of the Coulomb interaction (a second term in Eq. Q). In the 
case when spin moment operator S is well defined, so that density and potential matrices can be made diagonal in 
spin variables, it is convenient to define spin-polarized density and spin polarized potential. In a most general form 
exchange interaction allowing non-coUinearity of the spin polarization in space was introduced in density functional 
theory^ (DFT) by von Barth et alF'^^ and used later in Ref. I69ll70l : 

{[-V^ + Vextir) + / dVn(r')i rr + v^c{r\ n, rh)]ao + b^^{r- n, m) • S}-^,{r) = e,1',(r) (9) 

J \r -r' \ 

where 4'i(r) are one-electron spinors, Ui are Pauli matrices. Exchange-correlation potentials determined as variational 
derivatives over space charge density n(r) or magnetization vector m(r) are: 

5Exc[n,rh] - 5Exc[n,m] 5E.j,c[n,m] 

Vxc(r) = — — — — , h^c{r) = = m r). (10) 

on[r) om(r) dm 

The exchange interaction is described here by the local magnetic field b^c which acts on all orbitals in the same way. 
The potential matrix in LS basis for this field is diagonal in spin variables if quantization axis is chosen along direction 

of bxc- 

A strong spin-orbit coupling in 5/ shell produces an occupation matrix Eq. ^ and the corresponding potential 
matrix Eq. jnj which can not be made diagonal in spin variables by rotation of the quantization axis. Such rotation 
must be done to the direction of the spin moment vector, so that only Sz components were present while Sx,Sy 
components (responsible for off-diagonal in spin variable terms in Hamiltonian) set to zero. This is possible in the 
case where LS coupling scheme is valid and spin moment operator S is well defined, which is true for 3c? and 4/ 
elements compounds, but not for 5/ electrons of actinides. For later materials not LS coupling scheme is valid but 
jj coupling or intermediate coupling schemes where spin moment operator S is not well defined. 

That means that the widely used^S*^ form of spin-polarized potential based on Local Spin Density Approximation 
(LSDA) can not be applied for actinides, because LSDA neglects off-diagonal in spin variables terms in both spin 
density and potential. Those terms appear due to spin-orbit coupling Eq. ^ and suppress spin moment formation 
while exchange (Hund) part of Coulomb interaction prefers fully spin-polarized state with diagonal in spin variables 
occupation and potential matrices. A competition between spin-orbit coupling and exchange interaction determines 
the magnetic state of 5/ shell. Neglecting off-diagonal in spin variables potential terms in the methods using LSDA 
potential enhances the tendency to spin moment formation. (In Sec. IIIII we show that neglecting off-diagonal in 
spin variables potential terms in LDA-f-U-fSO method results in a very strong increase of calculated spin and orbital 
moment values.) This fact could explain why antiferromagnetic ground state with a large value of spin moment was 
obtained in fully relativistic spin-polarized method calculations-^'"*^ in disagreement with experimental data^*^*2i 



III. METALLIC PLUTONIUM 



A. 5 phase of plutonium 

Phase diagram of plutonium is the most complex of all elements,^-'- Of especial interest is fee 5 phase of Pu (a = 
4.636 A). A phase transition from S- phase to low-temperature monoclinic a-plutonium phase (a = 6.183 A, b = 
4.824 A, c = 10.973 A, and (3 = 101.79°) is accompanied with a very large (19%) volume contraction. In the present 
work we were interested in the problem of magnetic state for Pu metal in a and S phases. 

In the LDA-I-U calculation scheme the values of direct {U) and exchange {Jh) Coulomb parameters must be 
determined as a first step of calculation procedure. It can be done in ab-initio way via constrain LDA calculations^^*^ 
In our calculations Hund exchange parameter Jh was found to be Jh = 0.48 eV, much smaller than the value used in 
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FIG. 1: Partial f/^ and /^/^ 
LDA+U+SO method (b). 



contributions to 5/ band densities of states of (5-plutonium calculated with LDA (a) and 



previous LDA+U studies^ Coulomb repulsion parameter value U = 3.84 eV was obtained in a good agreement with 
previously used U = A eViS2*^2iSi22i2i24 Coulomb interaction parameters describe a screened Coulomb interaction 
between 5/ electrons and hence crucially depend on the channels of screening taken into account in constrain LDA 
calculations. That is true for the direct Coulomb parameter U, but exchange Coulomb parameter Jh corresponds 
to the difference of interaction energy for the electrons pairs with the opposite and the same spin directions. As the 
screening process is defined by the charge but not spin state of the ion, the screening contribution is cancelled for 
exchange Coulomb interaction and parameter Jh does not depend on the screening channels choice. For example for 
id elements compounds, direct Coulomb interaction parameter varies from 8 eV for late transition metal oxides (NiO 
and CuO) to 3 eV for the beginning of the row (Ti and V oxides). At the same time exchange Coulomb parameter 
Jh is in the range of 0.85-0.95 eV for the whole id row. In our constrain LDA calculations we have found that Jh 
value is equal 0.48 eV within the accuracy 0.01 eV for both metallic Pu phases and for all plutonium compounds 
investigated in the present work. 

In our constrain LDA calculations it was supposed that / shell is screened only by s, p, and d electrons, so that 5/ 
electrons were considered to be completely localized and not participating in the screening. LDA calculations with 
spin-orbit coupling taken into account show (Fig. ^a)) well separated subbands corresponding to the total moment 
values j = 5/2 and j = 7/2 with the Fermi level crossing the top of j — 5/2 band. This splitting of 5/ shell into 
nearly filled /^/^ and empty /^/^ subshells can justify an additional screening of /^/^ electrons by /"^^^ orbitals. In the 
constrain LDA calculations with this screening channel taken into account, Coulomb parameter value was obtained 
as [/ = 0.44 eV. 

Plutonium is considered to be on the border between completely localized 5/ electrons of americium and itinerant 
nature of them for early actinides, so that both limits, completely localized and totally itinerant, are not appropriate 
for Pu 5/ electrons. The correct value of U should be somewhere in between the values calculated in those two 
limits (3.84 eV without any / electrons participating in the screening and 0.44 eV for full /^/^ 

-screening). As an 

additional requirement to determine U value we have chosen an equality of the calculated equilibrium volume of 
fcc-Pu to the experimental value of 5 phase (see Fig.|2l. This requirement is fulfilled for U = 2.5 eV and we used this 
Coulomb interaction parameter value in LDA-I-U+SO calculations for metallic Pu in both a and 6 phases and for all 
Pu compounds investigated in this work. 

The LDA-I-U-I-SO calculations for metallic Pu in S phase gave a nonmagnetic ground state with zero values of spin 
S, orbital L, and total J moments. The occupation matrix Eq. (|2Jl has six eigenvalues close to unit (see in Table |2 
the row corresponding to Jh = 0.48 eV) and is nearly diagonal in the {jmj} basis of eigenfunctions of total moment 
J. That gives a configuration of Pu ion in jj coupling scheme. In the Fig. QJb) the partial densities of states for 
/^/^ and /^/^ subshell orbitals are presented. In agreement with the occupation matrix analysis one can see nearly 
completely filled /^/^ band with the Fermi level on the top of it and an empty /^^^ band. The separation between 
the centers of these bands is w 4 eV. 

The origin of this nonmagnetic LDA+U+SO solution can be traced to the results of standard LDA calculations 
without Coulomb correlation correction but with spin-orbit coupling taken into account (see Fig. ^a)). The 5/ band 
is split by spin-orbit coupling into well pronounced /^^^ and /^/^ bands with a separation between their centers w 
1.5 eV. Fermi level is in the pseudogap between subbands closer to the top of /^/^ band. Comparing Fig. ^a) and 
Fig. n^b) one can see that taking into account Coulomb correlations via LDA+U correction potential (Eq. O does 
not change qualitatively the band structure, only the separation between subbands increases from 1.5 eV to 4 eV 
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FIG. 2: Total energy as a function of volume for fcc-Pu calculated with various Coulomb parameter U values (Va corresponds 
to experimental volume value for Pu in 5 phase). Arrows indicate minima of the curves giving equilibrium volume values. 
Calculated equilibrium volume is equal to experimental value for the curve with U — 2.5 eV. 



TABLE I: Electronic configuration of Pu ion 5/ shell in (5-plutonium calculated with LDA+U+SO method {U — 2.5 eV) 
for various values of Hund exchange parameter Jh- The largest values of occupation matrices off-diagonal elements (OD) in 
{LS} and {jruj} basis sets are given in the second and third columns. The seven largest eigenvalues of occupation matrix are 
presented from fourth to tenth columns. Columns from eleven to thirteen show the calculated values for spin (S), orbital (L), 
and total (J) moments.-™" Last four columns contain partial contributions of and /® configurations and types of coupling for 
5/ shell of plutonium ion (see the text for explanations). 
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0.49 
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0.040 0.889 0.890 0.912 0.912 0.913 0.914 


0.083 


0.099 


0.016 


0.01 


0.99 


0.97 


0.03 


0.50 


0.433 


0.127 


0.041 0.880 0.884 0.909 0.918 0.918 0.925 


0.582 


0.692 


0.110 


0.04 0.96 


0.81 


0.19 


0.55 


0.412 


0.279 


0.045 0.838 0.884 0.912 0.929 0.938 0.945 


1.369 


1.640 


0.271 


0.11 


0.89 


0.54 


0.46 


0.60 


0.392 


0.346 


0.050 0.823 0.888 0.916 0.935 0.949 0.956 


1.746 


2.064 


0.318 


0.13 


0.87 


0.41 


0.59 



TABLE IL Electronic configuration of Pu ion in 5-plutonium calculated with LDA+U+SO method as a function of the 5/ 
band shift value A. (Notations are the same as for Table 0) 

A, eV OD{j^g} OD{j,„.} largest eigenvalues 5^ L J Jj LS 

0.040 0.889 0.889 0.912 0.912 0.912 0.912 1.00 1.00 0.00 

0.040 0.879 0.879 0.906 0.906 0.906 0.906 0.023 0.029 0.006 1.00 0.99 0.01 

0.040 0.864 0.866 0.897 0.900 0.903 0.905 0.210 0.286 0.076 0.03 0.97 0.93 0.07 

0.042 0.642 0.865 0.919 0.919 0.934 0.942 0.920 1.700 0.780 0.31 0.79 0.70 0.30 

0.043 0.254 0.910 0.924 0.942 0.954 0.959 1.316 3.040 1.724 0.69 0.31 0.56 0.44 

0.040 0.086 0.886 0.886 0.938 0.947 0.962 1.378 3.609 2.231 0.89 0.11 0.53 0.47 
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0.14 
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0.012 


0.27 
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0.044 


0.41 


0.427 


0.188 


0.82 


0.432 


0.256 


1.36 


0.419 


0.256 



according to the value of U — 2.5 eV. Another effect of Coulomb interaction is nearly pure orbital character of /^/^ 
and bands in LDA+U+SO calculations comparing with a significant admixture of /^/^ orbitals to the nominally 
/^/^ band and vice versa in the LDA results (Fig.nja)). One can say, that nonmagnetic Jh — solution with /^/^ 
subshell filled with six electrons and empty /^/^ band is already 'preformed' in LDA. Role of Coulomb correlations 
in LDA+U method is to make it more pronounced with a pure orbital nature of the bands and increased energy 
separation between them. 

Photoemission spectroscopy is a direct tool for investigating electronic structure. Total density of states obtained 
in our LDA+U+SO study were multiplied with the Fermi function corresponding to the temperature of experiment 
and broadened with 0.2 eV Lorentian to account for the instrumental resolution (see Fig. O. With the exception 



7 




Binding energy [eV] 



FIG. 3: Experimental (PES spectra for S- and a-plutonium with hv — 41 eV and T = 80 and angle-integrated PES for 
PuSb at 15 K— -) and calculated (thin line) photoemission spectra for metallic Pu in 5 and a phases and for PuSb. 

of the sharp peak near the Fermi level an overall agreement between experimental and calculated spectra for d-Pu 
is satisfactory. This peak is usually considered as a si gn of strong many-body effects, which can be described by 
Dynamical Mean-Field Theory53,54,55 (dmFT) (see Ref. l58l59() but not by static mean-field approximation which is a 
basis of LDA-I-U method. However, LDA+U can reproduce the lower and upper Hubbard bands, on which partially 
filled / band is split by Coulomb interaction. Indeed the peak in experimental spectrum at « 1 eV corresponding to 
the lower Hubbard band is well described by the calculated spectrum. 

Our results are in agreement with experimental data showing the absence of ordered or disordered magnetic moments 
for plutonium in both 5 and a phases^*^*^ However, all previous electronic structure calculations gave a magnetic 
(usually antiferromagnetic) state lower in energy than a nonmagnetic state. To clarify this problem we have carried out 
an investigation of the stability of our nonmagnetic ground state to the parameters of the calculations. The influence 
on it of the different approximations for exchange interaction and spin-orbit coupling terms in the Hamiltonian was 
investigated as well. 

As a first step we have varied the value of exchange Coulomb parameter Jh around Jh — 0.48 eV value obtained 
in our constrained LDA calculations. In Tabled the results of LDA+U+SO calculations with Jh values in the range 
from 0.40 eV till 0.60 eV are presented. One can see that even 15% increasing of Jh value is enough to result in a 
magnetic ground state with a sizable values of spin and orbital moments. However, the total moment value remains 
very small, less than 0.3 /is. 

In addition to moment values, we show in the Table values of the largest off-diagonal elements of occupation 
matrix Eq. ^ in {LS} and {jrrij} basis. As one can see, for the value Jh = 0.48 eV obtained in constrain LDA 
calculations, {jnij} basis is the most appropriate with largest off-diagonal elements values less than 0.01, while in 
{LS} basis these values are very large (more than 0.4). That can be interpreted as a pure case of jj coupling with a 
well defined total moment (J) but not spin (S) and orbital (L) moments. With increasing of Jh value the off-diagonal 
elements in {jnij} basis grow fast and for Jh = 0.60 eV become comparable with the off-diagonal elements in {LS} 
basis making jj coupling scheme not valid any more. As the off-diagonal elements in {LS} basis are still very large 
(w 0.4), the LS coupling scheme is also not appropriate and intermediate coupling scheme is needed to describe the 
magnetic solution. 

The eigenvalues of occupation matrix Eq. ||2J) presented in Tabled show six orbitals with occupancies close to unit 
for all values of Jh- Considering this one can conclude that Pu ion is predominantly in configuration in magnetic 
as well as in nonmagnetic state. Only nonmagnetic ground states corresponding to the values of Hund parameter 
Jh < 0.48 eV can be interpreted as a pure configuration in jj coupling scheme. Magnetic solutions for Jh > 
0.48 eV correspond to predominantly /® configuration in intermediate coupling scheme. 

In order to estimate, how close is a calculated state to one of the pure states (/^ or configurations, in jj or LS 
coupling schemes) we propose the following simple formula using calculated values of spin (S), orbital (L), and total 
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FIG. 4: (a) Band structure of S-Pu (the width of the Une shows the contribution of the 5/ states to the particular band) from 
LSDA+U calculations with Coulomb parameters values U — 4 eV and Jh = 0.85 eV. (b) The same as (a) but for Jh = 0.48 eV. 
(c) The same as (a) but by LDA+U+SO method (17 = 4 eV and Jh = 0.85 eV). (d) Same as (c) but for Jh = 0.48 eV (see 
text). 

TABLE III: Calculated values for spin, orbital and total moments from LDA+U+SO calculations with off-diagonal in spin 
variables elements in potential matrix (Eq. ^ set to zero. For comparison the results of LSDA and LSDA+U calculations are 
presented^ 



Method 




S 


L 


J 


LDA+U+SO {U = 


2.5 eV, Jh = 0.48 


eV) 2.43 


3.57 


1.14 


LSDA+U (ReL 33) 




2.50 


3.70 


1.20 


LSDA+U (ReL 59} 




2.55 


3.90 


1.45 


LDA+U+SO {U = 


Jh = 0.48 eV) 


1.70 


2.30 


0.60 


LSDA (ReL 4L) 




2.04 


1.75 


0.29 


LSDA (ReL 33) 




2.11 


1.94 


0.17 


LSDA (ReL^ 




2.25 


2.40 


0.15 



(J) moments. Total moment value is the same m both coupling schemes {jj or LS): J = for and J = 5/2 for 
If we have a mixed state {1 — x) ■ + x ■ then x can be defined as a; = J/2.5. Spin and orbital moment values 
for configuration are equal to: S" = 0, L = in jj coupling scheme and S' = 3, L = 3 in LS coupling scheme. 
For configuration S — 5/14 w 0.36, L — 20/7 « 2.86 in jj coupling scheme and S = 5/2, L = 5 in LS coupling 
scheme. One can define a mixed coupling scheme with a contribution of jj coupling equal to y and of LS coupling 
correspondingly to (1 — y). In this scheme the calculated values of orbital and spin moment will be equal to 

L^x- (2.86 • zj + 5 • (1 - y)) + (1 - x) • (0 • y + 3 • (1 - y)) (11) 
S = x- (0.36 • y + 2.5 • (1 - y)) + (1 - x) • (0 • y + 3 • (1 - y) (12) 

These formulas allow to determine the coefiicient y. 

The values of coefficients x and y calculated in this way are shown in Table ^ For values of Hund parameter 
Jh < 0.48 eV they correspond to the pure configuration in 100% jj coupling scheme. With increase of Jh value 
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the contribution of LS coupling increases and for Jh close to 0.60 eV both coupling scheme gave approximately 
equal contribution demonstrating a clear case of intermediate coupling scheme. However the configuration is still 
predominantly with not more than 10% admixture of 

Not only the strength of exchange interaction parameter Jh is crucial for the magnetism of the 5/ shell in metallic 
Pu. The Fermi level is on the top of /^/^ band (see Fig.QJa)) and a small shift of relative energy position of 5/ and 
other states can lead to redistribution of electrons between 5/ and spd-bands and hence change the / configuration 
and magnetic state of Pu ion. To investigate this effect we run the LDA+U+SO calculations with a constant positive 
potential A acting on 5/ electrons, which should result in a charge flow from 5/ bands to spd-bands (see Table EJ. 
Already relatively small value of the shift (A = 0.41 eV) is enough to produce a well pronounced magnetic ground 
state with a sizable admixture of configuration (around one third). For the shift value equal to 1.36 eV the Pu ion 
is in a pure configuration with a very large values of spin, orbital, and total moments. 

These constrain LDA+U+SO calculations demonstrate that a relative position of 5/ and other bands is of great 
importance for the resulting configuration and magnetic state of Pu ion. Later we will show that in all investigated in 
this work plutonium compounds we obtained a magnetic ground state in contrast to the nonmagnetic one for metallic 
Pu. 

In order to study the influence of the different approximations for exchange interaction term in the Hamiltonian we 
have performed at first LSDA+U (electron density and LDA potential were used in spin-polarized forms) calculations 
with the same Coulomb interaction parameters values as in Ref. [s^ (exchange parameter value corresponds to Jh — 
0.85 and U — 4 eV). In this calculation the strong magnetic state with a configuration close to was obtained 
in agreement with results of Ref. [s^. The 5/ bands dispersion is presented in Fig. Qfa) (the width of the band 
curve is proportional to the contribution of 5/ states to the band). It agrees well with the corresponding figure in 
Ref. Decreasing of the Coulomb exchange parameter value Jh from 0.85 eV to the value 0.48 eV (obtained in 
our constrain LDA calculations) gave still magnetic solution with spin S = 2.5 but the 5/ configuration is now more 
close to than to (the Fermi level (Fig.^b)) is now above the top of the band which was cut by the Fermi level 
in Fig. El^a)). As a next step we have used LDA+U+SO calculation scheme (jll8|l with non-spin-polarized electron 
density and LDA potential but keeping Hund parameter Jh — 0.85 e\42Si2S (see Fig. EIc)). In this calculation we 
have found magnetic ground state {S = 2.5, L — 2.7, and J = 0.2) with a configuration close to And only the 
calculation by LDA+U+SO with Coulomb exchange parameter value obtained in constrain LDA calculations {Jh = 
0.48 eV) gave (see Fig.0Jd)) a nonmagnetic state with S' = 0, L = 0, and J = 0. 

In the Sec.|n]we have discussed that the use of spin-polarized potential based on LSDA is equivalent to neglecting 
off-diagonal in spin variables exchange potential terms and hence enhances the tendency to spin moment formation. 
To check this we have performed LDA+U+SO calculations with the off-diagonal in spin variables terms in LDA+U 
potential (Eq. set to zero {U = 2.5 eV and Jh = 0.48 eV). As a result we indeed obtained a magnetic ground state 
(see Table llll|) which is close to the results of LSDA+U calculationsiSiS To compare with the results of fully relativistic 
spin-polarized calculation we have performed LDA+U+SO calculations with the off-diagonal in spin variables terms 
in LDA+U potential (Eq. O set to zero using the value of Coulomb parameter U = Jh = 0.48 eV. In Ref. it was 
shown that U — Jh can be regarded as an effective Coulomb interaction parameter Ueff in LSDA+U calculations, 
so the choice of [/ = Jh is equivalent to Ueff = 0. Such calculations gave solution with a large values of spin and 
orbital moments (see Table 1111}) similar to LSDA results. 

B. Q phase of plutonium 

A transition from 5 to a phase can be described as a volume contraction on 19% and a monoclinic distortion of the 
fee lattice. To separate the influence on the electronic structure and magnetic properties of the volume contraction and 
lattice distortion, we performed the LDA+U+SO calculation for fee lattice with the volume per Pu atom corresponding 
to a phase. The result is the same nonmagnetic solution: S ^ 0, L = 0, and J = in configuration as was obtained 
for i5 phase. A variation of Hund parameter Jh has shown that for Pu in fee-structure with the volume of a phase 
nonmagnetic solution is a little more stable than for 6 phase volume (while in the latter case already 15% increase of 
Jh was enough to produce the magnetic state, for a phase volume 20% increase is needed). 

LDA+U+SO calculation for real monoclinic crystal structure of a phase (see Fig. [SJ also gave a nonmagnetic 
ground state with S — 0, L ~ 0, J = in configuration. Comparing with the results for 5 phase (Fig. ^ one can 
see that bands become more broad (Fig.EIa)) due to the smaller volume and hence increased hybridization strength. 
Correspondingly /^/^ and /^/^ bands have a much more mixed character (a strong admixture of /^/^ orbitals to the 
nominally /^/^ band and vice versa) comparing with a more pure band character for S phase (Fig. ^a)). However, 
after taking into account Coulomb interaction in LDA+U+SO calculations (Fig. [S^b)) the band structure becomes 
very similar to the case of S phase (Fig. ^b)). /^/^ and /^/^ bands have now pure orbital character. /^/^ band is 
nearly completely filled with the Fermi level on the top of it and an /''/^ band is empty. The separation between the 
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FIG. 5: Partial f^'^ and 
LDA+U+SO method (b). 



/^/^ contributions to 5/ band densities of states of a plutonium calculated with LDA (a) and 



TABLE IV: Electronic configuration of Pu ion in various plutonium compounds calculated with LDA+U+SO method. Notations 
are the same as for Table in 



Compound ODf^gj 


0D{, 


PuN 


0.436 


0.239 


PuCoGas 


0.426 


0.226 


PuRh2 


0.348 


0.302 


PuSi2 


0.441 


0.142 


PuTe 


0.441 


0.253 


PuSb 


0.453 


0.320 



^.y largest eigenvalues 



L 



J 



r r 



LS 



0.118 0.162 0.923 0.930 0.960 0.965 0.978 1.279 3.339 2.060 0.82 0.18 0.57 0.43 
0.029 0.796 0.883 0.897 0.922 0.939 0.941 1.144 1.594 0.450 0.18 0.82 0.62 0.38 
0.127 0.525 0.610 0.883 0.908 0.924 0.927 1.546 3.421 1.875 0.75 0.25 0.46 0.54 
0.070 0.907 0.916 0.939 0.949 0.960 0.962 0.692 0.829 0.137 0.05 0.95 0.77 0.23 
0.026 0.895 0.904 0.908 0.928 0.940 0.946 0.907 1.069 0.162 0.06 0.94 0.70 0.30 
0.046 0.127 0.959 0.961 0.973 0.981 0.982 1.583 3.650 2.067 0.83 0.17 0.44 0.56 



centers of those bands is « 4 eV. 

The calculated photoemission spectra (see Fig.O is not very different from the corresponding curve for the S phase 
except that it is more smooth and does not show high and lower energy shoulders. The sharp peak near the Fermi 
level which is stronger in experimental spectrum for the a phase than for 5 phase spectrum is missing in the calculated 
spectrum as it was in the case for the S phase. Again we expect that Dynamical Mean-Field TheorjiSi^iS (DMFT) 
(see Ref. Is^lsgf) can solve this problem. But the peak in experimental spectrum at « 1 eV corresponding to the lower 
Hubbard band is well described by the calculated spectrum. 

As we have shown above, the calculations scheme presented here is different from the ones used in previous investi- 
gations of plutonium electronic structure and magnetic properties. A resulting nonmagnetic ground state (Pu ion in 
configuration with jj coupling scheme) was never obtained in any other calculations. In order to check validity of 
our method we have performed calculations for a series of Pu compounds. The main aim of this calculations was to 
investigate magnetic properties for these compounds and compare the results with the available experimental data. 



IV. PLUTONIUM COMPOUNDS 



PuN 



Plutonium mononitride crystallizes in a rocfc-saZt-type structure with a = 4.905 Ai2i Neutron diffraction showed 
no long-range order and magnetic moments larger than 0.25 From the magnetic susceptibility and specific heat 

measurements an antiferromagnetic transition at T^r = 13 K was proposed.''^ According to another magnetic suscep- 
tibility curve, PuN is a Curie- Weiss paramagnet with /ie// — 1.08 /is.^" Electronic structure of PuN was calculated in 
L(S)DA with and without relativistic correction!^ SIC-LSDA method showed valency +5 to be energetically preferred 
in the model with partially localized 5/ electrons of Pu ion^S 

Fig.EJa) andEJb) show total and partial / density of states of PuN obtained in our LDA with spin-orbit coupling 
calculation (when the Coulomb interaction correction was 'switched off'). Nitrogen 2p-states are strongly hybridized 
with Pu-/ states and the Fermi level is located inside the j — 5/2 subband. While j — 5/2 and j — 7/2 subbands 
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FIG. 6: (a) Total and partial densities of states of PuN calculated in LDA. (b) Partial f^^^ and /^''^ contributions in Pu 5/ 
band for PuN from LDA calculations, (c) The same as (a) calculated by LDA+U+SO. (d) The same as (b) calculated by 
LDA+U+SO. 



TABLE V: Values of effective paramagnetic moments (in Bohr's magnetons /is) calculated from the total moment value J 
- -■ — ' LS coupling schemes /if"// and l^^jj and their weighted value fiT/i in comparison with experimental data Me// (^^^ 



text). 



Compound J jj LS /ifff Me// Me// M^fT" 

PuN 2.060 0.57 0.43 2.16 0.72 1.54 1.08^ 

PuCoGas 0.450 0.62 0.38 0.69 0.23 0.52 0.68^ 

PuRha 1.875 0.46 0.54 1.99 0.66 1.27 0.882^ 

PuSi2 0.137 0.77 0.23 0.34 0.11 0.23 

PuTe 0.162 0.70 0.30 0.37 0.12 0.30 

PuSb 2.067 0.44 0.56 2.16 0.72 1.35 l.Oi 



are well separated from each other (see Elb)), they do not have pure orbital character with a strong admixture of 
j = 5/2 states to the formally j = 7/2 subband. The results of LDA+U+SO calculation are presented in Fig.EIe) 
and|Bfd). In contrast to metallic Pu case Coulomb interaction correction not only has increased energy separation 
between j = 5/2 and j — 7/2 subbands, but also led to the splitting of j = 5/2 subband into occupied and empty 
states (see the peak of j = 5/2 character just above the Fermi energy onE^d)). As this split-off band contain one 
electron per Pu ion, one can conclude that our results correspond to the configuration in agreement with formal 
Pu valency +3 in PuN. 

In Table Cxl (the first row corresponds to PuN) the calculated values for spin, orbital and total moments are 
presented. The total moment J = 2.218 is close to an ideal value j = 5/2 for configuration. Five largest eigenvalues 
of the occupation matrix also show a well defined configuration. The analysis of occupation matrices off-diagonal 
elements show that in {jrrij} basis the values of those elements are nearly two times smaller than the corresponding 
values for {LS} basis. That corresponds to the intermediate coupling scheme closer to jj coupling. The analysis 
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based on the calculated values for spin, orbital, and total moments (see Eq. (I12II in Sec. lIII^ gave 89% and 11% 
configuration with 57% of jj and 43% of LS couplings. 

In order to compare our results with the experimental magnetic measurements data we should obtain the magnetic 
moment value using calculated values of spin (S), orbital (L), and total moments (J). In the cases when spin-orbit 
coupling is weak and Russell-Saunders coupling {LS coupling) scheme is valid with S and L operators well defined, 
this is a very simple task: total magnetic moment value can be calculated as: 

Mtot = 2-S + L (13) 

However for strong spin-orbit coupling when or intermediate coupling schemes should be used, this problem becomes 
much more complicated. For the intermediate coupling scheme there is no general solution of this problem and for 
jj coupling it can be solved only for the free ions in pure configuration. An effective paramagnetic moment obtained 
from susceptibility measurements using Curie- Weiss law can be calculated as: 

Me// = g ■ ^/J■iJ + l)) ■ ^lB (14) 

The problem is to define Lande (7-factor which can be calculated for pure and configurations in LS or jj coupling 
schemes. As for configuration total moment value J = 0, one need to calculate y-factor for configuration only. 
For ground state of configuration in jj coupling scheme Lande factor gjj = 6/7 ~ 0.86. In LS coupling scheme 
its value is g^s = 2/7 « 0.29. As the latter value is nearly three times larger than the former, gjj and gi^g can give 
only an upper and lower limits of g-factor for the case of intermediate coupling. 

In Table W\ the values of effective paramagnetic moments calculated using Eq. H14() are presented. In addition to 
^J'^ff and Mej^ calculated using Lande factors gLS and gjj correspondingly, we have calculated their weighted value 
/ig^^'j using relative weight of LS and jj coupling obtained from Eq. H12|l . 

Results for PuN (the first row of the Table IVj) gave /^e/'/ — 1-^4 /zs in a reasonable agreement with the experimental 
value /^gfj — 1.08 /ib. 



B. PuCoGas 

Since first reported in 2002 by Sarrao et al. f^ superconductor PuCoGas has attracted a special attention being the 
first Pu-containing superconductor in well-known superconductor class of '115' materials Its superconducting 
temperature Tc = 18.5 K is an order of magnitude higher than for other '115' superconductorsi^^^ Specific heat 
coefficient value 7 = 77 mJ mol^^ K^^ above Tp and Curie- Weiss behavior of magnetic susceptibility with fi^ff — 
0.68 /iB are indicative of unconventional superconductivity of PuCoGas Tanaka et alM- used periodic 
Anderson model to describe d-wave superconducting state of PuCoGas. Its isostructural counterpart PuRhGas shows 
analogous properties but with the lower value of SC transition temperature Tc — 8 Ki2^ 

PuCoGas crystallizes in tetragonal PA/mmm space group with a = 4.232 A and c = 6.782 Electronic structure 
calculation of PuCoGas was performed in fully relativistic full-potential method of local orbitalsf^'^ In that work the 
paramagnetic state has a total energy value substantially higher than FM and AFM solutions. 

The LDA calculation without spin-orbit coupling®'^'^^ showed the Fermi level crossing Pu 5/ bands manifold without 
any splitting. Another calculation by relativistic LAPW method^l gave Fermi surface analogous to CeMIns series. 
Microscopic model for compounds with / ions in jj coupled state was extended to PuCoGas in R-of. 96. Antiferro- 
magnetic order was proposed in electronic structure calculation in Ref. l97l For our knowledge, in this work we report 
the first electronic structure calculation for PuCoGas with Coulomb interactions taken into account beyond LDA or 
GGA approximations. 

In PuN only 5/-bands are present at the Fermi level while fully occupied N-p band is situated substantially lower 
(Fig-EIa))- In contrast to that in PuCoGas Go-d and Ga-p bands cross the Fermi level in addition to 5/ bands of Pu 
(see Fig. [2{a)) so that 5/ states contribute only one third to the of density of states value at the Fermi energy. This 
fact leads to a very complicated general band structure for this compound (Fig.lHIa)) but the partial 5/ densities of 
state presented on Fig.Ob) shows well pronounced /^^^ and /^/^ subbands with « 1.5 eV separation between them. 
This picture is very close to the results for 5 phase of metallic Pu (Fig. H^b)) but with a larger overlapping of /^/^ 
and /^/^ subbands and the position of Fermi level a little deeper inside the /^^^ band. This is an effect of a stronger 
hybridization of 5/ orbitals with Co-d and Ga-p states in PuCoGas comparing with metallic Pu. 

The Coulomb interaction correction in LDA-I-U-I-SO calculations (Fig. (B^c) and (d)) results in a larger energy 
separation between /^/^ and /^/^ subbands (till 4 eV) and nearly pure orbital character of these bands, again very 
similar to metallic Pu (Fig. ^^d)). In Table HVl (the second row corresponds to PuCoGas) the calculated values for 
spin, orbital and total moments are presented. Pu ion for our solution is predominantly in configuration but with 
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FIG. 7: (a) Total and partial densities of states of PuCoGas calculated in LDA. (b) Partial f^^'^ and f'^^'^ contributions in Pu 
5/ band for PuCoGas from LDA calculations, (c) The same as (a) calculated by LDA+U+SO. (d) The same as (b) calculated 
by LDA+U+SO. 



a significant (18%) admixture of configuration. Sizable values of spin and orbital moment together with relatively 
small value of total moment gave similar to PuN case an intermediate coupling scheme closer to jj coupling. 

We would like to note, that in spite of very similar to metallic Pu band structure for 5/ states (compare Fig.[7Jb) and 
Fig. ^b)) in the LDA calculations, LDA+U+SO calculations scheme gave a magnetic ground state for PuCoGas in 
agreement with experimental datai^Si^Si^LS^iSS In Table El (second row) the values of effective paramagnetic moments 
calculated using Eq. (|14|1 are presented. The weighted value iJ-'^ff — 0.52 /is is smaller than experimental Me// — 

0.68 fiB, which is very close to the ^J■''Jff — 0.69 fiB calculated with Lande g-factor gjj obtained in jj coupling scheme. 



C. PuRha 

PuRh2 crystallizes in Laves structure (C15) a = 7.488 Ai2^ From the magnetic susceptibilitji^S PuRh2 is a Curie- 
Weiss paramagnet with peff = 0.88 m_b, 9p ~ -49 K, and Xmax = 4.4 x 10""^ emu mol~^. Resistivity agrees well 
with the magnetic susceptibility results.^* According to the specific heat results with 7 ~ 145 mJ mol~^ and 9d 
= 190 K, PuRh2 was classified as 'middle-weight' fermion system without temperature dependence of 7 in the region 
of low temperaturesjSS 

The dominant contribution to the PuRh2 band structure gives a broad partially filled Rh-d band (Fig. |SJa)). Pu 
5/ states (Fig.|SJb)) show significant hybridization with Rh-d band. However the general feature common to metallic 
Pu and all Pu compounds investigated in this work, a separation into /^/^ and /''/^ subbands still can be seen here 
(Fig. IHIb)). The position of Fermi level inside the /^/^ band similar to PuN case (Fig. Elb)) shows that resulting 
configuration should be close to Indeed LDA+U+SO calculations gave (Fig.EJc) and (d)) the /^/^ states split 
into occupied and empty bands. In contrast to PuN case (Fig. |SJd)) empty /^/^ band is not narrow but rather broad 
and has a two peak structure due to the strong hybridization with Rh-d states. 
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FIG. 8: (a) Total and partial densities of states of PuRli2 calculated in LDA. (b) Partial j^^"^ and /^/^ contributions in Pu 5/ 
band for PuRli2 from LDA calculations, (c) The same as (a) calculated by LDA+U+SO. (d) The same as (b) calculated by 
LDA+U+SO. 



In Table Hvl fthe third row corresponds to PuRh2) the calculated values for spin, orbital and total moments are 
presented. In agreement with the aforementioned density of states analysis, the Pu ion is predominantly in 
configuration with a significant admixture (25%) of configuration. A large values of spin moment results in 
intermediate coupling scheme closer to LS coupling. In Table (third row) the values of effective paramagnetic 
moments calculated using Eq. (|14|l are presented. The weighted value Mgj'/ — 1.27 /is is larger than experimental 
J = 0.88 [iB which still hes in the hmits of [x'^.r = 1.99 [Ib and = 0.66 \ib closer to LS couphng value [i^jj- 



D. PuSi2 

Structural and niagnetic properties of PuSi2 were reported by Boulet et alJ^ although this phase was earlier 
investigated in Ref. llOll This compound crystallizes in tetragonal ThSi2-type (space group lAi/amd) structure with 
a = 3.9707(3) A and c = 13.6809(5) Aiffi PuSi2 susceptibility curve shows Curie- Weiss behavior^°°i^°^ with p^fj = 
0.54 iiB, 0p = -58 K, and xo = 2.3 x 10"^ emu mol~^ fRef. llOOf) . Almost field independent resistivity shows a broad 
maximum at 18 K suggesting strong spin fluctuationsi^ffi 

The LDA band structure (Fig. El^a)) shows a broad deep pseudogap with a very small density of states value on 
the Fermi level. This pseudogap separates /^/^ and /^^^ subbands (Fig.l^b)) and position of the Fermi level exactly 
in the pseudogap gave a completely filled f^^^ band and hence a pure configuration could be expected. Indeed 
LDA+U+SO calculations gave (Fig.|^c) and (d)) a solution with the same type of band structure as without Coulomb 
interaction correction except increased energy separation between /^^^ and /^/^ subbands. 

In Table Hvl fthe fourth row corresponds to PuSi2) the calculated values for spin, orbital and total moments are 
presented. Pu ion is in nearly pure configuration with a small (5%) contribution of configuration. Relatively 
small values of spin and orbital moments gave an intermediate coupling scheme very close to pure jj coupling. The 
weighted value of effective paramagnetic moment /^e/'/ ~ ^.23 /is (Table 0) is smaller than experimental A*gj^ = 
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band for PuSi2 from LDA calculations, (c) The same as (a) calculated by LDA+U+SO. (d) The same as (b) calculated by 
LDA+U+SO. 
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E. PuTe 



PuTe crystallizes in NaCl-type structure with a — 6.183 A. Resistivity shows narrow-R ap se miconductor behavior 
and an electronic specific heat coefficient has a high value 7 = 60 mJ mol^^ K^^ (Ref. Il03j) . In another workiS^ 
the value of 7 = 30 mJ mol"^ was reported. While formal valency of Pu in PuTe is equal +2, an intermediate 
valent Pu+^-Pu+'^ state in this compound was proposed>i2S4i2Sii2Li2& Magnetic and optical properties of PuTe show a 
number of peculiarities >i22iii2iiii A structural phase transition from NaCl to CsCl phase was proposed from resistivity 
measurements 

Electronic structure of PuTe was calculated by Dirac equation-corrected ASAjiH relativistic LAPW,^^* and LSDA 
augmented pla ne w ave method (ASW).^^^ From these calculation, /^/^ and /^/^ subbands are split on 1 eV^^^ (or 
0.3 eV in Ref. Ill4l) . Results strongly depend on spin-orbit coupling strength.^^^ Based on the Mossbauer spectra 
resultSfii^ magnetic transition going from PuSb to PuTe with the vanishing of local moments was suggested. 

Results of calculations by LDA and LDA-I-U-I-SO are presented in Fig.llUland in Tablc llVl fthc fifth row corresponds 
to PuTe). The electronic structure and magnetic state are very similar to those for PuSi2, only the pseudogap is 
not so well developed. Nearly pure configuration corresponds to the formal Pu valency +2 with a small values of 
magnetic moments described by 30% admixture of LS coupling to jj coupling scheme. The arguments against valency 
+Z m Ref. Eos 106 107 loj were based on the very large ionic radius value for Pu+^ ion. However, in spite of the 
formal configuration in our results, the total number of 5/ electrons is equal to 5.65 (effect of strong hybridization 
of 5/orbitals with Sb-p states), which corresponds to the intermediate valency Pu+^-Pu+'^. Another interesting fact 
is that going from PuTe to PuSb the number of 5/ electrons decreases from 5.65 to 5.16, two times smaller than 
expected from formal valency difference. A very small value of magnetic moments obtained in our calculations (Table 
IV|) agrees well with nonmagnetic state of Pu ions in PuTeiii^ 
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FIG. 10: (a) Total and partial densities of states of PuTe calculated in LDA. (b) Partial /^''^ and /^'''^ contributions in Pu 
5f band for PuTe from LDA calculations, (c) The same as (a) calculated by LDA+U+SO. (d) The same as (b) calculated by 
LDA+U+SO. 



F. PuSb 



From magnetization measurements and neutron experiments an antiferromagnetic spin structure was found with 
a strong <001> anisotropy and ord ering temperature 75 KiiiLii^ Preferable intermediate coupling type closer to 
iS'-typc was proposed in Ref. '119* (model of wave function consisting of 90% LS- and 10 % jj-coupled wave 
functionsyi. gave a possibility to AFM transition leaving strong anisotropy) . Neutron-diffraction experiment reveals 
low-temperature Pu ion magnetic moment = 0.75 /is perpendicular to ferromagnetic (001) planes with Neel temper- 
ature Tjv = 85 K and first-order transition to the incommensurate ferromagnetic phaseii22iii2ii2£ Effective magnetic 
moment of paramagnetic phase from susceptibility was estimated as /i = 1.0 Model of localized 5/ electrons 

with weak hybridization to conducting electron states was proposed»i2Sii2S Electrical resistivity shows metallic be- 
havior with Kondo-like broad maximum at 106 K.^^^ X-ray photoelectron spectroscopy (XPS) and high-resolution 
valence-band ultraviolet photoelectron spectroscopy (UPS) showed strong localized nature of 5/ Pu electrons in PuSb 
and pronounced configuration.^^"' 

Results of calculations by LDA and LDA-I-U-I-SO are presented in Fig.lllland in Table Hvl ft he sixth row corresponds 
to PuSb). The electronic structure and magnetic state are very similar to PuN, according to the formal valency +2 
for Pu in both compounds. 

We have calculated photoemission spectrum for PuSb (see Fig. I^J. An agreement with experimental cmvei^ is a 
good one. The main peak at « 1.5 cV corresponding to Pu-5/ states (see Fig. Illf cll is reproduced quite well. 

In Table IVlf sixth row corresponds to PuSb) the values of effective paramagnetic moments calculated using Eq. H14|) 
are presented. The calculated fi^ff = 1.35 is close to experimental value fxlYf = 1-00 Ms- 
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FIG. 11: (a) Total and partial densities of states of PuSb calculated in LDA. (b) Partial /^^^ and /^/^ contributions in Pu 
5/ band for PuSb from LDA calculations, (c) The same as (a) calculated by LDA+U+SO. (d) The same as (b) calculated by 
LDA+U+SO. 



V. CONCLUSION 



In this work we present the results of theoretical investigation for Pu ion magnetic state in metallic plutonium and 
plutonium compounds. In contrast to all previous theoretical studies but in agreement with experimental measure- 
ments we have found for metallic Pu in both a and S phases a nonmagnetic ground state with configuration of 5/ 
shell in pure jj coupling scheme. A strong spin-orbit coupling for 5/ electrons results in a splitting of / bands into 
well separated /^/^ and /^/^ subbands with /^/^ band nearly fully occupied and /^/^ band empty giving 'preformed' 
configuration. Taking into account Coulomb interaction via LDA+U potential does not change this nonmagnetic 
ground state. 

We have shown that approximately equal strength of spin-orbit coupling and exchange interaction, whose matrices 
can not be made simultaneously diagonal in the same basis set, does not allow the use of the simplified diagonal forms 
of the corresponding Hamiltonian terms. That is true also for spin-polarized potential used in LSDA. Only a general 
non-diagonal matrix form of exchange interaction is appropriate for 5/ electrons. 

We have calculated also a series of plutonium compounds with different formal valency and the calculated magnetic 
moments values agree reasonably well with experimental data. 

Comparison of the calculated and experimental photoemission spectra show that while LDA+U+SO method based 
on static mean-field approximation can reproduce Hubbard bands appearing in such strongly correlated materials 
as plutonium and its compounds, the quasiparticle peak on the Fermi energy can not be described by this method. 
The more elaborate Dynamical Mean-Field Theory (DMFT) is needed here. The results obtained in LDA+U+SO 
calculations could be used a basis for further DMFT studies. 
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